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The dynamics of molecular motors which occasionally detach from a heterogeneous track like 
DNA or RNA is considered. Motivated by recent single molecule experiments, we study a simple 
model for a motor moving along a disordered track using chemical energy while an external force 
opposes its motion. The motors also have finite processivity, i.e., they can leave the track with a 
position dependent rate. We show that the response of the system to disorder in the hopping off rate 
depends on the value of the external force. For most values of the external force, strong disorder 
causes the motors which survive for long times on the track to be localized at preferred positions. 
However, near the stall force, localization occurs for any amount of disorder. Existence of localized 
states near the top of the band implies a stretched exponential contribution to the decay of the 
survival probability. To obtain these results, we study the complex eigenvalue spectrum of the time 
evolution operator. A similar spectral analysis also provides a very efflcient method for studying 
the dynamics of motors with infinite processivity. 



I. INTRODUCTION 



Single molecule experiments which study molecular motors provide a powerful tool for understanding their function 
0. By applying a mechanical force, one can often discern details of their reaction steps. For example, a recent paper 
considered the motion of a kincsin moving along a microtubule under the influence of a force opposing its motion 
0. By measuring the relationship between the applied force, F, and the velocity of the kincsin, v, at different ATP 
concentrations it was possible to infer that the chemical reaction cycle contains at least one load dependent transition. 

More recently, it has beenpossible to study, using similar techniques, the motion of RNA polymerase (RNAp) 
DNA polymerase (DNAp) helicases |^ and A-exonuclease • In contrast to kincsin and myosin, which move along 
homogeneous polymer filaments, these motors move along DNA, which is inherently a heterogeneous track, with the 
energy landscape determined by the nucleotide sequence |^. Indeed, the dynamics of these motors seems far richer 
than that of kinesin Q or, say, myosin V which move along periodic filaments. Theoretically, it is expected that 
heterogeneous and homogeneous (or periodic) tracks can give rise to very different motor dynamics. On long time 
and large length scales the motion of motors moving along a homogeneous track is described well by a random walker 
moving along a tilted potential. In contrast, the motion of molecular motors which use chemical energy to move along 
disordered filaments is described by a random walker moving on a random forcing energy landscape [Tol llll | . The 
effective energy difference between two points separated by m nucleotides scales as y/m. The large energy barriers 
implied by such landscapes lead to anomalous dynamics when the overall tilt of the energy landscape is small; 
the displacement of the motor grows as a sublinear power of time. For molecular motors this corresponds to applying 
an external force strong enough to place the system near the stall point of the motor. 

The theories described above do not treat the effect on the dynamics of a (position-dependent) detachment proba- 
bility of motors from the track. We will refer to models where no unbinding of the motor from the track is allowed as 
"infinite processivity" models. In this paper we study theoretically the effect of a finite processivity on the behavior of 
molecular motors. Previous work has considered the effect of detachment mainly on the dynamics of many molecular 
motors or with homogeneous tracks jl,'^ll4ilT5lllfiLll7j . We consider a motor which can leave the polymer on which it is 
moving but can never rebind to it. This case is relevant to single molecule experiments which follow a specific motor or 
a dilute concentration of motors on a long DNA track. In such experiments (especially if a background hydrodynamic 
flow washes away detached motors) the probability of rebinding to the track after detachment is negligible. 

We show that when molecular motors are moving along a homogeneous (or periodic) filament the long time dynamics 
of those motors that remain attached is unaltered by detachment events. However, for motors moving along a 
heterogeneous track this is not the case. For proteins such as RNAp and DNAp that walk along DNA, the detachment 
rate depends on the monomer on which the motor is located. The disorder in the detachment rates is thus correlated 
with the disorder in the hopping rates as both are determined by the same DNA sequence. When the detachment 
rate varies from monomer to monomer, the motors' dynamics can be influenced by spatially localized eigenfunctions 
of the evolution operator. For the regime where the displacement of the motor is linear as a function of time, strong 
disorder causes the last motors which remain attached to the track to stall out before falling off. Moreover, when the 
displacement of the motor as a function of time is sublinear, any amount of disorder in the hopping off rates causes the 
last motors which remain attached to the track to halt before falling off. A schematic representation of the resulting 
"phase diagram" for the dynamics is shown in Fig. ^ 
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FIG. 1: Schematic behavior of the influence of disorder in the rates for leaving the track on the dynamics of the motors. A 
represents the strength of the disorder, measured through the variance of the hopping off rates divided by the square of their 
mean. Note that when A = the motors are always delocalized except at Fstaii. This is emphasized by the shading. It should 
be stressed that the notions "localized" and "delocalized" here refer to motors which remain on the track for long times. 



To show these results we use methods developed in [I8|,ll9j, in the context of the physics of vortices, and numerically 
study the eigenvalue spectrum of the non-hermitian evolution operator for the probability distribution of the motor. 
As will be shown, the characteristics of this operator's spectrum (which include localized and delocalized states as 
well as a mobility edge) allow the long time, large length-scale dynamics to be obtained in a straightforward manner. 
Interestingly, the dynamics of motors on heterogeneous tracks with finite processivity is similar to the non-hcrmitian 
statistical mechanics arising in models of vortex physics |l8l Il9j . 

To support these results further we also study analytically a toy model which consists of a directed walker among 
traps with a broad distribution of release times. The model was studied for the case of infinite processivity in We 
show that in the infinite processivity case the model yields the same spectrum as observed in the numerics. Moreover, 
we show that for a motor with finite processivity, when the corresponding infinite processivity model shows sublinear 
drift, any amount of disorder in the hopping off rates leads to localization of the cigcnfunction, consistent with our 
numerical results. 

The paper is organized as follows: In Section II the model we study is introduced and known results for the infinite 
processivity limit are reviewed. In Section III the infinite processivity model is analyzed using the eigenvalue spectrum 
of the evolution operator. It is shown that the long time, large length-scale behavior of the model can be extracted 
from the properties of the spectrum. In Section IV we study the spectral properties of an evolution operator that 
takes account of detachment of motors from the track and their physical implications. Specifically, the implications 
for measurements of the probability of finding a motor on the track as a function of time are discussed. Finally, in 
Section V the analysis of the toy model is presented. 



II. THE MODEL 

The model we consider was introduced and studied in the limit of infinite processivity in It is inspired by 

previous models of molecular motors but is simple enough so that exact solutions can be found, in 

the limit of infinite processivity, with and without disorder. The model is defined on a discrete lattice, x = 0, 1, 2 . . ., 
with distinct a (even) and b (odd) sites and a distance ag between lattice points. A monomer of the track (a nucleotide, 
say) is taken to be of size 2ao. The arrangement is shown schematically in Figs. [21 and 01 To model the two internal 
states of the model we take even sites to have an energy £ = while odd sites have an energy e ~ Ae. The transition 
rates depicted in Fig. Otake the form 

u;- = (ae^^'/^+c.)e-^^/^-^/2^ 
= [a + ijj)ef /"^^ 

= (a'e'^''/^ + ^')e-^^/^+^/'^ (1) 
= (a'+^')e^//2^, 

where we have set the Boltzmann constant to be Ub = 1- Note that there are two parallel channels for the transitions. 
The first, represented by contributions containing a and a', arise from utilization of the chemical energy A/i. The 
chemical energy difference could be, for example, a result of an excess concentration of NTP's (or just ATP) with 
respect to its thermal equilibrium value. The second channel, represented by the terms containing u and cj', correspond 
to thermal transitions unassisted by chemical energy. In addition, the externally applied force F , with / = i^ap, biases 
the motion. To model finite processivity we add rates w^jj and lu^g corresponding to the detachment from the track 
at even and odd sites respectively. 
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FIG. 2: Setup modelled. The motor is moving from the — end to the + end, driven by hydrolysis of nucleotide tri-phosphate. 
A force is pulling on the motor in the opposite direction. The track shown is made up of two types of monomers (depicted as 
shaded and white areas). Although this schematic suggests a microtubule with disorder in the protein constituents, we actually 
have in mind motors moving on DNA or RNA templates, with four distinct nucleotides. Elaborations which make our model 
more realistic at a microscopic level should not affect the predictions for long time, large length-scale dynamics. 
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FIG. 3: Graphical representation of the model for molecular motors. The distinct even and odd sites are denoted by a and h 
respectively. 



Before turning to the study of finite processivity we review some of the results which are known for the infinite 
processivity limit We show later that these results may be reproduced by studying the spectral properties of 

the evolution operator of the probability distribution. It is the spectral method that will allow us to determine 
the dynamics most readily when the processivity is finite. In the homogeneous case (a, a',w and lo' independent of 
position), when no disorder is present the model is described on long times and large length-scales by a random walker 
moving along a potential with an overall tilt between two even sites which is given by 

Ai; = rin 



' ^^""y («eAM/T + ^) (c,' + ,,') j + ■ 

We refer to such an energy landscape as an effective energy landscape. It is an alternative description of the dynamics 
associated with the rates of Eq. which in general will not satisfy detailed balance since they describe nonequilibrium 
processes. Note that when f — and the chemical potential difference A/i = 0, one has AE = and no net motion 
is generated. Also, when there is directional symmetry in the transition rates a ~ a' , co = uj' (reflecting directional 
symmetry in the DNA track) and / = one has AE = even when A/i ^ 0. Absent this symmetry, chemical energy 
can be converted to motion and an effective tilted potential is generated. Similar conditions for biased motion have 
been shown to exist for continuum models [22, . The effect of the externally applied force is simply to change the 
overall tilt in the potential. Thus, for motors moving along a homogeneous (or periodic) polymer, one expects the 
velocity to change continuously as / is changed (see Fig. ^J. 

For motors moving on heterogeneous filaments the situation is very different. Here the set of parameters {p} — 
{a, a', w, w', Ae} is drawn from a random distribution. Each set of parameters describes the dynamics on a given 
type of monomer. Using the results presented above it is easy to see using equation ||2Jl that the total effective energy 
change after m monomers is given by 



£;(m) = ^Ai?(0. (3) 

1=1 

Here, each AE{m) corresponds to an independent set of values of {p] drawn randomly for the m*'* monomer. Assuming 
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that AE(rn) is drawn from a random distribution with a finite variance, the effective energy landscape corresponds 
to a biased random walk. Such energy landscapes are typically referred to as random forcing energy landscapes. 

The above scenario applies as long as the chemical potential difference A/i 7^ 0. In the case when A/x = it is easy 
to see that E{m) = unless we allow for the energy at even sites also to vary and take the value £{m) |2^. In this 
case we obtain 



E{m) = 2fm + e{m) , (4) 

corresponding to a random energy landscape provided e(m) has only short range correlations. 

The dynamical behaviors of random walkers in random forcing or random energy landscapes have been studied in 
detail 0, . Using the results of Dcrrida [2^ one can calculate the transition points between the different regimes 
including the effect of randomness for our model jlll |. Upon denoting averages over the disorder rates by an overbar 
one finds the following regimes. 
Regime I: Ordinary biased diffusion occurs when 



/<-^ln(:U^j , (5) 
4 V wrw^ / 



4 V wtrwt 



where the subscript / = denotes that / has been set to zero when evaluating the average. In this regime (x) — vt 
and (x^) — (x)^ = 2Dt for long times, where the angular brackets denote an average over different thermal histories 
of a particle starting from a particular point. Another way of stating this result is that for large times the calculated 
velocity and diffusion constant, defined by the above relations, do not depend on the size of the time window, tyy, 
over which they are evaluated. This results only holds for biases satisfying Eqs. (O and (jH)). 

Regime II: The calculated diffusion constant now depends on the size of the time window, tw, over which it is 
evaluated. The velocity does not. Namely, in this region (x) = vt and {x^) — (x)'^ ^ t^/^, where 1 < < 2. In the 
infinite tw limit the diffusion constant, D = limt„,-^oo((2;(tn/)^) ~ {^itw))'^)/tw^ diverges. This anomaly occurs in 
the ranges 



T , / wtrwtr \ ^ „ . T 



In ^-4- </<_ , (7) 



and 



?lnP^ </<?lnp^y . (8) 

Regime III: Here both the velocity and the diffusion constant arc functions of the size of the time window tw, and 
{x) ^ t^ and (x^) — {x)^ ^ t'^^, where fi{f) < 1. In the infinite tw limit the velocity v = \imt„->oo{{x{tw) — x{0))) /tw 
vanishes. The diffusion constant, D = \hnt^^,^cc{{x{tw)'^) — {x{tw))'^)/tw, either vanishes or diverges depending on 
whether fi < 1/2 or fi> 1/2. This behavior occurs when 



_Z,J^!/r^] <f<liJ^^] . (9) 

Sinai diffusion: Exactly at the stall force fs, {x) — and (x^) ^ (ln(i/r))'*, where r is the microscopic time needed 
to move across one monomer. The "stall force" corresponding to a disordered track is defined by 



^.m ( ) . (10) 



2 V wirwt 



The resulting behavior as the force is varied is summarized qualitatively in Fig. ^in the limit of an infinite averaging 
time window tw ■ Most notable is the region of forces (region III) over which the displacement of the motor is sublincar 
and the usual long time velocity vanishes. Experiments are performed with finite tw- The measured velocity then 
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FIG. 4: Schematic behavior of the velocity for a heterogeneous linear motor track as a function of the applied force. It is 
assumed that chemical forces (from the NTP hydrolysis) lead to a positive velocity in the absence of a force. The anomalous 
dynamics arises in the vicinity of a stall force defined by Eq. 1101 The different dynamical regimes defined in the text are 
denoted in the figure. The striped line on the /-axis denotes the region where anomalous diffusion is present. 
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FIG. 5: The velocity as a function of f /T for different values of tw- Here Afi/T = 3 and parameters were chosen with 
equal probability to be either {p} — {5, 1, 0.3, 1, 0} or {p} — {4, 0.1, 0.7, 1, 0} (see text for notation). The calculated regime of 
anomalous velocity is 0.5116 < f /T < 0.699. Data incorporates 100 runs (thus averaging over thermal fiuctuations) for a single 
realization of the disorder. 



behaves in region III as t^^ ^, smoothing the curve shown in Fig. ^ and naturally giving rise to a convex shape of 
the velocity-force curve. This convexity is demonstrated in Fig. where the model was simulated and the velocity 
measured using different averaging time windows tw (see Appendix A for details). As is evident in this figure, the 
larger tw, the closer is the velocity- force curve to that shown in Fig. ^ In Fig. typical trajectories of the motor 
on the track are shown for different values of f /T for a single realization of the disorder. Plateaus and jumps appear 
as one moves closer to the region of anomalous displacement, and the dynamics is controlled by deep minima in the 
effective energy landscape with rapid transitions between them. Such motion is typical of random forcing energy 
landscapes 



III. INFINITE PROCESSIVITY: A SPECTRAL ANALYSIS APPROACH 



The behavior of the infinite processivity limit, as summarized above, is rather well understood. In this section 
we show that the dynamics with infinite processivity may also be deduced by an alternative approach: the spectral 
properties of the evolution operator of the probability density are exploited to deduce the long time, large length-scale 
properties of the model. The important features of the eigenvalues and eigenfunctions which characterize the long 
time dynamics are expected to be insensitive to details of the model. The method has previously been applied in 
the study of the physics of vortex lines in superconductors 0, ^ and population dynamics [23 . Earlier studies 
considered the spectral properties of a random walker on a random forcing energy landscape after "gauging away" 
the external bias, making the evolution operator hcrmitian [T^ . However, as we will show, a more direct approach 
reveals important features (such as complex eigenvalue spectra) which are not easily observed in the approach of ^3 • 
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Although here the approach is used to reproduce known results, it might prove useful in the study of more complicated 
models with an infinite processivity where an analytic solution is not possible. In addition, the spectral approach to 
the infinite processivity model will be very useful as a benchmark once we address the question of finite processivity. 

In what follows, wc will often consider a coarse-grained effective dynamics in which each new lattice site represents 
a unit cell, containing one a and one h site, of the original lattice (Eq. ^and Fig. OJ. The lattice spacing in this 
new coarse-grained model is a = 2ao. Below we first review the spectrum (and the associate dynamics) of the trivial 
homogeneous model before turning to the disordered one. 



A. Homogeneous Model 

As discussed above the homogeneous problem is well described on long time-scales and large length-scales by 
a random walker moving along a tilted potential. To see this from the model's spectrum, we first solve for the 
probability of being in odd sites and substitute the solution in the equation for the probability of being in even sites. 
The resulting coarse-grained equation |lll |. in the long-time limit, for the probability, P{x,t), of being at site x at 
time t is given by 

{w^ + w^)dtP{x,t) = w';;;*w^P{x -2,t) + w'^w^P{x + 2,t) 

~ {wT^r +w^w^)p{x,t) . (11) 

In the continuum limit, which describes the long time, large length-scale behavior of the model we have 

dtP{x, t) = DdlPix, t) - vd^Pix, t) , (12) 

with D ~ a^{w^w^ + w^w^)/2{w'j^ + w^) and v ~ a{w'j^w'^ — w^w^)/ (w^ + w^). Here, as noted above, 
a = 2ao is a lattice constant analogous to a unit cell size in conventional solid state physics. We consider a system of 
size L = Na and perform a Laplace transform in time and a Fourier transform in space so that 



P{x,t) = ^.g(/c)e''^^+^('^)* 
fc 



(13) 
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The coefficient g{k) will depend on the initial conditions. If periodic boundary condition are imposed on a lattice 
of size L, the wavevectors k specifying the (delocalized) eigenfunctions ^'^(a;) ~ e**^^ are quantized, fc„ = 27rn/L, 
n — 0,±1, .... In general, since the evolution operator is non-hermitian eigenvalues may be complex. The complex 
eigenvalues must come in complex conjugate pairs to ensure a real probability density P{x,t) [l9j. The eigenvalue 
spectrum is given by 

A(fc) = -Dk^ - ivk , (14) 

corresponding to diffusion with drift. The dependence of the real part of the spectrum on k describes diffusion, while 
the dependence of the imaginary part describes drift with a constant velocity. This identification is possible when 
periodic boundary conditions allow a current to flow in the system even in the limit oi t —^ oo. 

B. Heterogeneous Model 

For the heterogeneous model the situation is more subtle, because the Fourier transform cannot be used to diag- 
onalize the evolution operator. Nevertheless, we will show that a numerical analysis of the eigenvalue spectrum can 
be used to deduce the dynamical properties of the model. The role of the wave numbers in the homogeneous model 
is taken by winding numbers |28l |. 

The winding number, defined only for complex eigenfunctions, is given by the number of times the eigenfunction 
spirals around the origin in the complex plane as it traverses the whole lattice with periodic boundary conditions. For 
homogeneous systems and eigenfunctions which behave like e'*"'"^ we have /c„ = 2TTn/Na, implying a winding number 
n. For disordered systems, it can be shown [2^ that the winding number increases linearly with the eigenvalue index n 
which orders the eigenfunctions from the lowest value of |A(n)| to the highest. That is. the winding number associated 
with an eigenvalue of magnitude |A(n)| is ±n. Moreover, if the eigenvalue X{n) is in the upper part of the complex 
A plane its winding number is n. Its complex conjugate pair A*(n) will then have a winding number —n. For real 
eigenvalues the eigenfunctions can be chosen to be real and the winding number is not defined. In this case, the 
eigenfunctions can be classified by the number of zeros, as in one dimensional quantum mechanics |29l |. 

The eigenvalue spectrum of the heterogeneous model with periodic boundary conditions can be used to encapsulate 
the long time dynamical properties of the model. Specifically, the dependence of the real and imaginary parts of the 
eigenvalues on the winding number signals the properties of the mean square displacement and drift of the motor 
molecule within our model. 

To this end we consider the model of Sec. II with, as before, two types of monomers drawn at random. The 
set of parameters representing each type of monomer {p} = {a, a' ,lu,uj' , Ae} is chosen with probability 1/2 to be 
{-Pi} = {o^ii c«'ij i^i) "^1: ^£1} and with probability 1/2 to be {^2} = {0:2, 0-2^^2, ^'2 ^ ^^2}- The chemical potential 
difference A/i is assumed to be the same for both types of monomers. Allowing Afi to depend on the type of monomer 
does not alter the qualitative long time behavior of the model. 

We calculated the eigenvalue spectrum of the evolution operator numerically by diagonalizing a matrix with 
a specific realization of the disorder (see Appendix B for more details). From the exact solution of the model 
we expect three different regimes. Below a typical spectrum from each of the regimes is shown and examined 
in detail. Throughout we set T = 1 and use the parameters {pi} = {ai, a'j^, wi, , Aei} = {6,1,1,6,0} and 
{P2} = {^2, 0:2, W2, ^£2} = {1-2, 1, 1, 1-2, 0} (again with equal probabihty) and e^^/-^ = 10. The results are 
unchanged for similar sets of parameters. For these values of the parameters the exact locations of the transitions 
between the different regimes can be readily calculated using the results summarized in Sec. II. 

Regime I: For / < 0.2256 and / > 0.518 the motion is biased diffusion. A typical spectrum of eigenvalues in the 
complex plane (for / = 0) is shown in Fig. |7| We show only the eigenvalues with the smallest |A|, because those 
with larger values are non-universal and are highly dependent on the details of the model used to describe the motion 
of the motor. The spectrum has two branches in the complex A plane. Also, it has no purely real parts (with the 
exception of A = 0). 

For usual diffusion with uniform drift in a homogeneous system, as discussed above, A(fc) = —Dk^ — ivk. The 
imaginary part of the spectrum corresponds to drift, the real part to diffusion. A similar behavior might be expected 
in our heterogeneous model for A(fc) where k is related to the winding number as k = kn = 2-Kn/Na. This can be 
seen in a plot of the real and imaginary parts of the spectrum as a function of the eigenvalue index as defined above. 
Since in this regime the dynamics is biased diffusion we expect Im(A(n)) oc n, while Re(A(n)) cx . 

To test this hypothesis, we order the eigenvalues according to their magnitude and plot Im(A(n)) and Re(A(n)). 
Fig. IHlshows Im(A(n)) for small n. One can see that indeed the slope in linear. Fig. |2|shows — Re(A(n)) for small n 
on a log-log scale along with a line — Re(A(n)) = Av?' with A some constant. At small n indeed Re{\{n)) oc v? for 
more than a decade. 
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FIG. 7: The eigenvalues obtained for / = 0. Here the system size is A'' = 4500. Shown are the 140 eigenvalues with the lowest 
value of |A|. In the region near the origin we expect that Im(A) oc ±y^Ite(A)|. 




FIG. 8: The imaginary part of the eigenvalue, Im(A(n)) for / = as a function of n. Here the system size is A'^ = 4500. Shown 
are the 140 eigenvalues with the lowest value of |A|. 




FIG. 9: The real part of the eigenvalue, — Re(A) for / = as a function of the winding number n. Here the system size is 
A'' = 4500. Shown are the 140 eigenvalues with the lowest value of |A|. The solid line is the function An^ with A a constant. 
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FIG. 10: The eigenvalues obtained for / — 0.25. Here the system size is A'^ = 4500. Shown are the 140 eigenvalues with the 
lowest value of |A|. Near the origin we expect that lm(A) oc ±|Re(A)|^'''', with 1 < < 2. 
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FIG. 11: The imaginary part of the eigenvalue for / — 0.25 as a function of the winding number n. Here the system size is 
= 4500. Shown are the 140 eigenvalues with the lowest value of |A|. 



Regime II: Next we turn to look at the region where the diffusion constant depends on the size of the time averaging 
window used to evaluate it, while the velocity does not. For the parameters used in the numerics this occurs for force 
windows given by 0.2256 < / < 0.2882 and 0.4555 < / < 0.518. A typical spectrum of eigenvalues in the complex 
plane is shown in Fig. 1101 Again we concentrate on the small |A| part of the spectrum and examine Im(A(n)) and 
Re{X{n)). Fig. ^Jprcscnts a plot of Im(A(n)) for small n. The expected dependence on n is linear. A careful analysis 
of the curve shows that Im(A(n)) can be fit well to An + Bit? for small ?7 (here we have used that fact that Im(A(n)) 
is expected to be an odd function of n) . The coefficient B decreases as the size of the system studied is increased. 
Such a correction is expected if the time to reach the asymptotic behavior (x) ^ vt is so large that it is comparable to 
the relaxation time of the system. For larger systems the correction will become less important and the asymptotic 
behavior will be observed more easily. Fig. E| shows —Re{\{n)) for small n on a log-log scale along with the lines 
An and Bn^ with A and B some constants. At small n, i?e(A(n)) indeed does not behave as n^. It shows a smaller 
slope in accordance with a behavior consistent with 1 < /x < 2. In fact, from naive dimensional considerations (A 
behaves as \/t while x ^ n) this is the behavior expected when {x^) — (x)^ ~ t^/^. 

Regime III: Finally, we turn to look at the region where the velocity is effectively zero as — * oo and the diffusion 
is again anomalous. For the parameters used in our numerical analysis this regime occurs for 0.2882 < / < 0.518. 
A spectrum of eigenvalues in the complex plane is shown in Fig. 1131 Again we concentrate on the small |A| part of 
the spectrum. To get a clear picture one needs to plot both Im(A(n)) and Re(A(n)). The positive part Im(A(7i)) is 
plotted on a log-log plot in Fig. ^1 Here the displacement of the particle is expected to behave as with ^ < 1 . This 
leads to an expected dependence Im(A(n)) cx n^^^ . A reference line An with A a constant is plotted for comparison. 
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FIG. 13: The eigenvalue spectrum obtained for / — 0.31, where both diffusion and drift are anomalous. Here the system size 
is A'^ = 4500. Shown are the 40 eigenvalues with the lowest value of |A|. Near the origin, we expect Im(A) oc ±lRe(A)l. 



Clearly, Im(A(ri)) behaves as expected. 

Fig. 1151 shows — Re(A(n)) for small n on a log-log scale along with a line An^ with A some constant. As can 
be seen at small n indeed Re(A(n)) does not behave as n^. Here a similar argument as before leads one to expect 
— Re(A(n)) cx n^^^ for small n, consistent with the numerics. 

We comment that the analysis is only possible not too far from the transition point into regime 111. Deep inside the 
region with anomalous velocity, the small |A| spectrum becomes very noisy due to finite size effects. Analysis of the 
spectrum then becomes difficult due to these strong sample-to-sample fluctuations. Hence we have not attempted an 
analysis in the Sinai diffusion regime, where the energy landscape is nearly horizontal. 



IV. FINITE PROCESSIVITY 



A. Homogeneous Model 



As discussed above, the homogeneous model is well described on long time scales and large length scales by a 
random walker moving along a tilted potential. The unbinding of motors from the track adds a non-conserving term 
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FIG. 14: The imaginary part of the eigenvalue for / — 0.31 as a function of the winding number n. Here the system size is 
A*' = 4500. Shown are the 40 eigenvalues with the lowest value of |A|. A fit to Im(A) oc —\n\^^^, over only one decade, gives 
= 0.6 ± 0.1. The function An with A a constant is plotted for reference. 




FIG. 15: The real part of the eigenvalue, —Re{X) for / — 0.31 as a function of n. Here the system size is N = 4500. Shown are 
the 40 eigenvalues with the lowest value of |A|. A fit to Re(A) oc n^^'^, over only one decade, gives /i = 0.6 ± 0.1. The function 
An with A a constant is plotted for reference. 



to the equation for the probability density. In a continuum description, one now has 

dtP{x, t) = DdlP{x, t) - vd^P{x, t) - WoffP{x, t) . (15) 

In terms of the microscopic model as before one has D = a?'{w'j^w'^ + w^w^)/2{wi^ + w^) and v ~ a{w^w^ — 
wtw'a )/{w-^ + wD- In addition, here WoS = {w'^^siw^ + ^r) + <ff(w'r + wr))/iw^ + wf ). 

Again, we consider the model with periodic boundary conditions and perform a Laplace transform in time and a 
Fourier transform in space. The eigenvalue spectrum is given by X{k) = ^Dk- — ivk — WoS- When Woff — the 
motion is diffusion with drift. However, when Woff > the probability density decays exponentially to the empty 
track state P{x,t) = with a typical time scale 1/woff. Note, however, that the eigenfunctions are identical to the 
infinite processivity case. In fact, the probability density of the motor can be written as 

P{x, t) = e-"'°«-*P^„„=o(a;, t) . (16) 

Up to a rescaling of the probability density the time evolution is unaltered. That is, the motors which remain on the 
track are unaffected by a constant hopping off rate. 
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B. Heterogeneous Model 

If the hopping off rates are uniform along the track, but the remaining parameters are chosen randomly, w^g and 
w^g are the same for all types of monomers, and it is straightforward to see that the effect is the same as for a 
homogeneous system: the dynamics of the motors which remain on the track are unaffected by the non-conservation. 

We now show that, when randomness in the rates for hopping off the track is introduced into the system, the 
dynamics of the motors which remain on the track can be altered in a dramatic way. The eigenfunctions of the 
evolution operator are affected by the heterogeneous hopping off rates as well as by the sequence heterogeneity as it 
affects the local diffusion constant and drift velocity. As we shall see, disorder in the hopping off rates has a profound 
effect on the behavior of motors which remain on the track. We stress that, with applications to molecular motors in 
mind, the disorder in the detachment rates we consider is correlated with the local hopping rates. 

As for the conserving heterogeneous model with its random force landscape, we study the dynamics by considering 
the spectral properties of the evolution operator of the model. The effect of detachment is included by setting the 
rates w^g and w^g, corresponding to leaving the track at even and odd sites respectively, to be non-zero and site 
dependent. We find two types of behavior which depend on the drift properties of the motor for the corresponding 
conserving model. When the drift is linear in time the last motors to be left on the track localize only beyond a 
critical disorder strength in the hopping off rates. Similar behavior arises in the physics of vortices [TsL IT^ and in 
population dynamics |27l |. The behavior when the displacement of the motor for the conserving model is sublinear 
(i.e. when it is not possible to define a drift velocity) is different. Now, the last motors to remain on the track are 
always localized for any strength of the disorder. 

1. Ballistic displacement 

First, we consider the regime where the displacement of the motor in the conserving model is linear as a function of 
time {fi > 1). We will show that when the disorder in the rates of leaving the track is strong enough the eigenvalues 
with the smallest |A| become purely real. This is a regime of localized eigenfunctions and corresponds to motors with 
zero drift velocity. This feature is present over the whole region where fj, > 1. 

To explore this behavior, we have studied the eigenvalue spectrum of the model numerically for systems of size 
4500. In Fig. ^| results are shown for the case / = (/x > 2). However, the localization effects which are of interest 
to us here remain unchanged for all values of / such that fi > 1. The parameters used are the same as used previously 
in the paper augmented by w^g = 0.01 and w^g = 0.01 {w^g = 0.02 and w^g = 0.02) for the parameter set {pi} 
({^2})- As can be seen, except for the lowest |A| eigenvalue, all eigenvalues have an imaginary part. Thus, the motors 
are biased to move across the lattice in a positive direction. In Fig. El we show the spectrum with w^g = 0.01 and 
w^g = 0.01 {w^g — 0.08 and w^g = 0.08) for the set {pi} {{P2}), which corresponds to a larger disorder strength. 
Here eigenvalues beyond a certain mobility edge become real, implying that the long lasting motors are described by 
localized states [il,|2g. Note, that the disorder in the hopping off rates is correlated with the disorder in the hopping 
rates themselves, as required if both are due to the same underlying heterogeneous polynucleotide sequence. 

The transition between moving and localized motors can be understood, in the case of ballistic displacement, by 
adapting the ideas of Ref. [s^- We consider a simplified model for the motors. The energy landscape consists of two 
slopes si and S2 (corresponding to two types of monomers). The velocity on a track with a slope si (52) is denoted 
be vi {v2)- The rate for leaving the track on a region with a slope si is assumed to be w, while the rate for leaving 
the track on a region with slope S2 is 0. Note that as before the local detachment rates are correlated with the local 
biases. 

Next, compare two "survival" strategies for the motor in the long time limit. In the first it spreads across the lattice 
probing regions with slope si and S2- In the other it remains in regions with slope §2 where detachment rates are 
zero. The probability of staying in such regions can be estimated by solving for a random walker moving on a tilted 
potential of some size with absorbing boundary conditions at both ends. One finds |30| | that this probability behaves 
as 

Pstay(t) ^ e-''^'*/^^ , (17) 

where D is the diffusion coefficient on a track with slope S2. The probability of survival of a particle moving across 
the lattice can also be evaluated. One expects 

Pmovc(t) - e--^' , (18) 

where c is a constant which will depend on the disorder averaged velocity. For large enough w, the particles that 
choose the localized strategy clearly will have a better chance of survival. One therefore expects that as the disorder 
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FIG. 16: The eigenvalue spectrum for / = and /i > 2 with rates for hopping off the track as specified in the text. Except 
for the largest eigenvalue, all eigenvalues have an imaginary part. All eigenfunctions, moreover describe extended states with 
a well-defined winding number. Shown are the 140 eigenvalues with the lowest value of |A| and the system size is N = 4500. 
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FIG. 17: The eigenvalue spectrum for f — with rates for hopping off the track as specified in the text. As can be seen the 
eigenvalues with the lowest value of |A| become localized. Shown are the 300 eigenvalues with the lowest value of |At and the 
system size is A'^ = 4500. 

in the hopping off rates is increased a transition from localized states to delocalized states will occur at the high end 
of the spectrum (corresponding to small k and small winding numbers) . Note that the argument would not be altered 
if the rate for hopping off from negative slope regions were non-zero, although, the transition point between the two 
regimes would be shifted. 

Finally wc comment that this argument could also be used for a model with a random energy landscape (sec Eq. 
0}. This case is very similar to problems which have been studied in the context of the vortex physics |l8l Il9l| and 
population dynamics |26l| . For completeness. Appendix C presents eigenvalue spectra, for weak and strong disorder, 
in that case. 



2. Suhlinear displacement 

A typical spectrum in the regime of sublinear displacement, /.t < 1, is shown in Fig. 1181 Note the small band of 
localized states with real eigenvalues near the top of the spectrum (small |A|). Again the parameters used are the 
same as used previously with the faUing from the track rates w'^^ = 0.01 and w^jj = 0.01 (w^fj = 0.02 and w^fj = 0.02) 
for the set {pi} {{p2])- Note, that even for these small values of the falling off rates localized states appear. This 
effect persists even for falling off rates smaller by more than an order of magnitude and for various values of / in 
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FIG. 18: The eigenvalue spectrum for / = 0.31 with rates for hopping off the track as specified in the text. The smaU band of 
real eigenvalues at the top of the spectrum characterizes localized states. Shown are the 140 eigenvalues with the lowest value 
of |A|, and the system size is again A*' = 4500. 
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FIG. 19: The eigenvalue spectrum for / — 0.31 with rates for hopping off the track as specified in the text. The spectrum is 
qualitatively similar to that of Fig. I18l with weak disorder. The band of real eigenvalues at the top of the spectrum characterizes 
localized states and now contains, as expected, more states than Fig. 1181 Shown are the 140 eigenvalues with the lowest value 
of |A|, and the system size is again A'' = 4500. 



this regime as long as large enough systems were studied. For comparison Fig. 1191 shows the spectrum in the case of 
strong disorder. Here w^g = 0.01 and w^g = 0.01 {w^g = 0.06 and w^g = 0.06) for the parameter set {pi} {{p2})- As 
always in this paper, the two parameter sets arc assumed t ooccur with equal probabihty. As can be seen from the 
figure, while details of the spectrum arc different from the weak disorder case shown in Fig. 1181 thcv are qualitatively 
the same. 



3. Implications for the probability of finding a motor on the track as a function of time 

Next, we discuss the effect of the localized states occurring near the top (low |A|) of the spectrum on the probability 
of finding a motor on the track as a function of time. It may be possible to reveal these localized states by looking 
at the density decay as a function of time of a dilute concentration of fluoresccntly labeled motors placed on a track. 
The probability of finding a motor on the track as a function of time, Ps{t), could then be monitored by looking at 
the decay of the fluorescent signal as a function of time (possibly averaging over several experiments). 
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The implication of localized states for the probability Pg (t) has already been considered in the context of random 
walkers subject to the influence of random traps |30l |. The arguments are unaltered for the case of a random forcing 
energy landscape, and here we consider a simple version. A more detailed proof can be carried out along the lines 
of Ref. [25. We consider a system where the hopping and detachment rates can only assume two values. We 
are interested in the behavior in the long-time limit. When localized states exist the eigenfunctions are strongly 
peaked in regions where the probability of detachment from the track is small. The probability for such a region of 
length I to occur behaves as e~^'' where 7 is a constant. On the other hand the time to leave the region behaves 
as g-'"2*/'^D+xDt/i ^ -^^jjere x is a constant of order unity and we have included higher order corrections to Eq. El 
associated with the finite size I of the region. Note that we have used the fact that the detachment rate is correlated 
with the local hopping rate so that the behavior inside the region is diffusive irrespective of the value of /i. Summing 
over the contribution from regions of different lengths one obtains 

Here io depends on the constants 7 and x and D the diffusion coefficient. As expected the leading behavior is still 
exponential. However, a signature of the localized states appears in the correction which has a stretched exponential 
form. 

In contrast, when no localized states exist the behavior is dominated by the detachment rates. In that case no 
diffusive corrections to the detachment rate from a region of size I are present. The stretched exponential correction 
to the decay of displayed in Eq. ((T^ will be absent. 

In summary, in the region of sublinear displacement one expects a stretched exponential behavior to always be 
present at long times when the probability of finding a motor on the track is monitored. In the ballistic displacement 
regime the correction will be present only for strong enough disorder in the detachment rates. 



V. TOY MODEL 



The understanding of the dynamics of random walkers (with infinite proccssivit y) o n random forcing energy land- 
scapes has been enhanced by a simple toy model introduced by Bouchaud et. al. |l2l |. The model builds on the fact 
that at finite tilt the sojourn time r at any site is found to have a broad distribution ^(r) ~ r"^^"*"''-* for large t. This 
behavior suggests that the dynamics could be mimicked by a directed walk between traps with a broad release time 
distribution. Specifically, the model consists of a particle moving unidirectionally on a lattice with hopping rates Wk 
between site k — 1 and k which are drawn from a probability distribution which satisfies \I'(W^) w AW^~^ as ^ 0. 
Although the equivalence to the original model cannot be justified rigorously (since back stepping is ignored), the two 
models are known to exhibit the same long time behavior. 

Here we analyze the model from a somewhat different perspective and show that indeed it yields for the infinite 
processivity limit the types of spectra described in Section HI. Moreover, we study the finite processivity model in the 
limit of weak disorder in the rates for leaving the track. Here the rate for leaving the track at site fc, Wk, is assumed to 
be uncorrelated with the hopping rate Wk- If one views the directed model as a coarse-grained version of the original 
undirected model, then one expects that the hopping rates Wk arc determined by the depth of the traps in the energy 
landscape. In contrast, since the particles spend most of their time near the bottom of the traps, the coarse-grained 
off rates in the directed model will depend primarily on the local off rates near the trap bottoms in the original model, 
suggesting that the Wk and Wk should indeed be uncorrelated. With this simplification, we show that similar to the 
more complex models discussed earlier in this paper, any strength of disorder in the rates for leaving the lattice causes 
the eigenfunctions describing the long-time behavior of the model to be localized for ji < 1. In contrast for /i > 1, 
weak disorder leaves the eigenfunctions delocalized. 

It is convenient to first consider a more general model with non-zero rates for leaving the track. We derive an 
equation for the eigenvalue spectrum of the model and then analyze the spectra separately in the infinite and finite 
processivity limits. 

The time evolution of the model is described by the set of master equations 

= WkPk-i - Wk+iPk - Wk+iPk ■ (20) 

dt 

After a Laplace transform, the corresponding equations for the eigenfunctions is then 



(A + Wk+i + Wk+i)Pk{X) = WkPk-i{X) , 



(21) 
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which yields a recursion relation for the Pk{X) which can be readily solved. The possible eigenvalues A are then 
determined by imposing periodic boundary conditions, yielding 

1 . n 



where the last line holds in the limit when the system size N is taken to infinity and the overline denotes as before 
average over disorder. Care must be taken when A is real. In this case the recursion relation might allow for some 
A + Wk+i + Wk to be zero. In that case Pk{X) can assume any value (to be set by normalization), while Pk-i{X) = 0. 
Such eigenfunction will be referred to as "localized" , with a detailed justification given below. When A is complex, 
we refer to the eigenfunctions as delocalized. 
According to Eq. we need to calculate 



Since we are considering the long time behavior and small rates for leaving the track we consider the limit jA + wj ^ 1. 
It will be useful to hold w fixed for now and average over it later. We denote the function obtained by averaging only 
over W in Eq. H23|l by Q{X, w). We find that it depends on the value of the exponent /x as follows: 
Case 1: ^ > 2. Here one finds to leading order in (A + w) that 



Case II: 1 < /i < 2. Here W^'^ diverges and one has 



Q{X,w) = ~(^-^yX + w) + B{X + wr , (25) 

where B = A7r//i sin(7r/i). Here, A is the amplitude of the tail of the probability distribution, ^'(ly) « AW^^^ . 
Case III: ^ < 1. Here both W~'^ and W~'^ diverge and one has 

Q{X,w) = B{X + wY ■ (26) 

With these results at hand we now turn to analyze the infinite proccssivity limit. 



A. Infinite Processivity limit 

In this case we set w = for all sites. To analyze the eigenvalue spectrum we rewrite Eq. H22|l in the form 

R{X) = 2'Kin/N . (27) 

Consider the solution for the different cases 
Case I: /i > 2. Here one has 



2WA^ — (^)a+(^)^. (28) 

It is straight forward to solve the quadratic equation for A and realize that Im(A) cx n while Re(A) (x in agreement 

with the numerics of section IhH 

Case II: 1 < /i < 2. In this regime one has to solve 



2^T^n/N = -(^^Yx) + B{Xr , (29) 

where B ~ A7r//isin(7r/i). Again, a straightforward analysis shows that Im(A) oc n while Re(A) oc in agreement 
with the numerics. 

Case III: /i < 1. Finally, in this case one has 

2Tiin/N = B{XY , (30) 
which evidently gives Im(A) cx -n}/^ and Re(A) cx n^/^. 
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B. Finite Processivity 

Next, we turn to the finite processivity case. Here it will be sufficient to consider only the leading order contribution 
to R{X). To simplify, we study the case where w can take only two values wi — w and 1^2 — 0, which occur with 
probability pi and p2. The general case can be analyzed very similarly. This assumption gives, after performing the 
average over w, 

R{E)^PiQ{E,w)+p2Q{E,0) . (31) 

As in the previous subsection we analyze behavior of the solutions of the equations for A and w small. Here it will 
only be necessary to consider the cases /i > 1 and fj, < 1. 
Case I: /i > 1. Here we have 

~{^){pi{X + w)+p2X)^27Tni/N (32) 

with n an integer. Clearly, to solve the equation A must have an imaginary part. (Note that we are missing the 
dependence of the real part of the solution on n since we neglected higher order terms in A). An imaginary part of 
the eigenvalue implies that delocalized states exist for small disorder in w |l9j |. 
Case IT. fi < 1. Here one has 

B{p,{X + wr+P2Xn = ^ (33) 

Consider first the case n = 0. In this case the eigenvalue must be real [sj and since the decay of particles from the 
system can not be faster than e~™* one must have —w < A < 0. Clearly, such a solution can not exist. Thus, the 
assumption leading to Eq. H33|) . that none of the Pk(X) is zero, fails. As we argued early, Pfe(A) = away from a 
particular site, can only occur when A is real. It is possible to show that in this case the density is peaked near the 
site where Pk{X) = decaying exponentially fast as k increases. We conclude that for any strength of disorder in the 
hopping off rates the eigenfunctions become localized. The lack of solution for cases with n > can be proved by 
expanding around the lowest |A| (real) eigenvalue and discovering that a solution is impossible. 
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APPENDIX A: SIMULATIONS 



In the following wc describe briefly the procedure we used to simulate the model Eq. To make the simulation 
efficient we first normalize the rates so that the largest one is equal one. Then, at each step we choose with equal 
probability between moving the motor to the right or left on the lattice. Following this choice a random number 
is drawn from a uniform distribution. The motor is moved in the chosen direction provided the random number is 
smaller than the corresponding rate. This protocol ensures relaxation to equilibrium in the absence of chemical or 
mechanics driving forces. Time in the simulations is measured by the number of attempted moves. 



APPENDIX B: EVOLUTION OPERATOR 



Here we describe the construction of the evolution operator corresponding to model Eq. First consider the 

case when the rates for leaving the track are set to zero. To do this we write the Master equation describing the time 
evolution as 

dt\P{t)) = M\P{t)) , (Bl) 

where \P{t)) is a vector with components Pn{t), the probability of being at site n at time t |3l|. Ai is the evolution 
operator of the model with non-zero components given by 

Mi4+i = Mi^i-i = for i even 

Mi^i+i = Mi^i-i = w7 for * odd 

Mi^r = -Mi-i,^ - Mi+i,t . (B2) 
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FIG. 20: The eigenvalue spectrum for the random energy model with rates for hopping off the track as specified in the text. 
Here an example of the weak disorder regime spectrum is shown. Note that the right most part of the Re(A) axis is displaced 
slightly below zero. Shown are the 140 eigenvalues with the lowest value of |A| and the system size is A'^ = 4500. 



The last relation ensures conservation of probability. For a heterogeneous model the rates corresponding to a given 
type of monomer A^i.i+i, A4i+ij, Mi+2,i+i, M.i+i,i+2 with i even are chosen at random. The diagonal terms are 
then automatically given by Eq. Ilj2l To study the eigenvalue spectrum periodic boundary conditions are imposed. 
Finally, the hopping off rate from site i, hi (given by w'^^ or w^^j depending on wether i is even or odd) is added 
through Aii^i A^ j,i — hi. The spectra presented in this paper were calculated using MatLab. 



APPENDIX C: RANDOM ENERGY MODEL 



In this Appendix we discuss the spectrum of the random energy version of the model, described by Eq. ^ which also 
allows the energy at even sites to vary. The rates are slightly modified so that they to satisfy local detailed balance 
with respect to these energies. An example of the relevance of random energy models to biological system can be found 
in [33, which studies the diffusion of regulatory proteins along DNA. Another realization (although there is no analogy 
of "falling off" ) is translocation of ssDNA through a pore with identical environments on the trans and cis sides [ll| ■ 
As stated in the text, the long time lar ge w avelength properties of this model are very similar to those of a model 
studied in the context of vortex physics We use model (QJ with T = 1, A/i = 0, T = 1 and use the parameters 
{pi} = {ai, a'l, wi, a;i,ei, e'l} = {6, 1, 1, 6,0, 1.5} and {^2} = {^2, ^2, ^3, £2, £2} = {1-2, 1, 1, 1-2, 0.5, 0.8}. Here Ci, 
e'i denotes that energy at even and odd sites respectively. 

In Fig. 1201 and I2n the spectra of the evolution operator of the model are shown for weak (w^g = 0.1 and w^jj = 0.1 
«ff = 0.12 and wl^ = 0.12) for the parameter set {pi} {{P2})) and strong disorder (w^fj = 0.1 and = 0.1 
(w^g = 0.4 and w^^j = 0.4) for the set {pi} ({^2}) in the detachment rates. As can be seen from the figures strong 
enough disorder in the hopping off rates, as before, causes eigenvalue with small |A| to be real and therefore associated 
with localized eigenfunctions [l9l |. 
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FIG. 21: The eigenvalue spectrum for the random energy model with rates for hopping off the track as specified in the text. 
Here an example of the strong disorder regime spectrum is shown. Note that the right most part of the Re(A) axis is displaced 
slightly below zero. Shown are the 300 eigenvalues with the lowest value of |A| and the system size is A'' = 4500. 
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